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ABSTRACT 


The object of this research was to model an actual 
structure, Bullard Hall, in order to devise a more energy 
efficient means to control the internal temperature of that 
Structure. The research included a study of all major 
thermodynamic forces acting on the structure, including 
solar radiation. Once the model was developed it was com- 


pared with the actual structure to determine its usefulness. 
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t. ST ROoUCT ION 


To design and implement an optimal control scheme, the 
designer must first understand the dynamics of the system to 
be controlled. Usually the control engineer works with the 
transfer function of the system; either derived from the 
mipelse response or by computing the transfer function from 
a model of the system. In the study of thermal structural 
analysis, it is extremely difficult to devise a simple 
transfer function due to the complexity of the system and 
Ene Verrbesey of the forces acting on a structure. To obtain 
the impulse response of an actual structure would require 
the expenditure of large amounts of energy and the develop- 
ment of some set of standards with which to judge the 
response, this method of obtaining the transfer function was 
Semsuadered, but proved too difficult to obtain useful data. 
With the advent of the computer, the engineer can now 
COnsenuce models of the system using equations to describe 
the system and its environment. These equations can later 
be used to determine which parameters or forces have the 
most effect on the system, and by considering only the major 
forces acting on the system a smaller model and eventually a 


usable transfer function can be derived. 





The objective of this research was to develop a working 
model of the thermodynamic system consisting of the internal 
environment of Bullard Hall and the forces acting on that 
Sumleoure which would affect the internal environment. This 
Sumuevure Was selected for its simplicity and availability, 
and once the building was modeled, the model could be used 
to devise more cost effective means to conserve energy while 
maintaining the comfort of the personnel in the structure. 
Due to the complexity of the modeling process, the time 
involved, and the lack of sufficient recorded data an actual 
working control scheme was not developed. 

In Chapter II the appropriate energy equations for the 
Structure are stated and simplified to make the problem 
realizable. All the major sources of energy loss and gain 
are discussed and modeled for inclusion in an overall system 
model. Particular emphasis is made to model the solar 
energy portion of the system in an attempt to better utilize 
that particular energy resource. 

Ciapees £11 1S a compilation of the formulas derived in 
Chapter II and a computer run using the model compared with 
dewleleddta. Chapver Ili also explores the use of cost 
functions and how they would be used to design an optimal 
control system to reduce the use of nonrenewable energy 


sources while keeping the building comfortable. 





Chapter IV discusses tne problems encountered in the 
collection of usable data for the model and the instrumenta- 
tion used. This chapter also provides suggestions for 
future work in the area of instrumentation development and 
control studies in the area of energy conservation. 

The final chapter summarizes the results of this 
investigation and presents conclusions and suggestions that 


might be considered for further study. 





II. THEORY AND THE MODEL 


The» object of a model is to portray the appearance, 
nature, or performance of a device or system. It is assumed 
that with sufficient information a set of mathematical 
equations or formulas can be derived which will allow the 
modeler to emulate any particular component or collection of 
components of a system. Once a model has been created and 
verified to actually emulate the real system, the modeler 
can then change parameters in the model and study the effects 
of these changes. If a particular response is desired from 
the actual system, the model can be used to test the effects 
of changes on the system before costly changes are made. 

The use of a model to study tne effects of changes in a sys- 
tem is part of the science of controls, and if a specific 
effect or overall condition is specified, with some neces- 
sary tradeoffs, the process then becomes one of finding the 
optimum control. This model is used to describe the thermo- 
Gvaameec ertects of the building (Bullard Hall) by the air, 
internal equipment, and solar radiation surrounding or 
acting on the structure. Once verified, a suitable cost 
function or set of values will be selected to find the 
optimum control for the building's internal temperature. 


The major concern is to obtain maximum comfort for the 





people who must occupy the building while using the minimum 
amount of energy. In order to simulate the thermodynamic 
properties of a building, certain information must be 
collected, simplified, and condensed into a model. The 
first problem is to define the system and its boundaries. 
For the purposes of this model the building was considered 
to be a constant volume system. This assumption leads 
immediately to the energy conservation equation for a 


constant volume [1], the equation is: 


Jeo Mie + ee = M(m + ee + AE/dt 


W. = Rate of Shaft Work or energy transferred to 


the system as work in BTU/hr. 


Q = Rate of energy transferred to the system as 
edu a uy ar. 
M = Mass flow across the Svovcenmooundary adn. lor ar. 
(e+Pv) = Energy contained in the mass crossing the 
system boundary in BTU/1b. 
AE/dt = Change in energy storage term in BTU/hr. 


At this point some simplifying assumptions are made: 

1) There is no energy transferred to the building as 
work done on the system boundary. 

2) The PV term for both the incoming flow and exiting 


flow is equal. This assumption is reasonable, since it 





requires approximately the same amount of energy to get into 
the building as it does out. 

3) The energy that is stored in the building is stored 
only as internal energy; there are no mechanical or chemical 
SvOrage units in the building. 

4) The energy storage medium of the building is its 
mass, not the people, machinery, or furniture in the 
building. This last assumption allows us to express the 


energy storage term as: 
meee = Mo xeCp x di/dt 


mieGesim 1S the total mass of the building, Cp is the 
Specific Heat capacity of the mass substance and dT/dt is 
the rate of temperature change of the building's mass. 


The original equation now reduces to: 


Q + M(e,, -e P= Mx Cp ey de 


out 


At this point some further simplifying assumptions are 
made. 

5) The energy storage medium of Bullard Hall will be 
the concrete mass of the building. 

6) The air inside the building will be assumed to be in 
Cie~mabeequilibrium with the internal concrete walls. This 
assumption allows one to take air temperature in the build- 


ing as if it were the temperature of the mass. During data 


10 





acquisition this assumption proved to be quite accurate for 
the air temperature for particular regions of the building 
differed from the wall temperature by no more than 2.0 
degrees. 

7) Since the system is so massive, it possesses a very 
Slow response time, and therefore data were taken hourly. 


Heat loss and gain {Q + M(e, - Cate) } of any system can 


n 
occur in three ways: conduction, convection, and radiation. 
These three methods will be discussed later in the form of 
emvelope loss, infiltration, and solar gain. It is impor- 
tant to note that all the heat loss and gain terms are 
considered together. During the 1960's and 70's, most 
literature and textbooks on building heating and cooling 
only considered solar heat gain as a load on the air- 
Segarcloning system and did not consider it as a contributor 
in heating a building. In the particular case of Bullard 
Hall, the solar gain is substantial and can not be ignored. 
The actual mass of the building is the sum of the 
machinery, air, and the physical structure. In this model 
we will only consider the mass of the structure since the 
total mass of the air plus the machinery in the building 
together amounts to less than 2% of the building mass. For 
a low mass structure the air and machinery may be a signi- 
RuccwiemOuemvity, DUG in this case it is not. the building 


maSS waS calculated using the blueprints [2] and the 1952 


a 








ns 





Dui idingespeciiications [3]. The following table is a 


listing of those calculations: 


ITEM NUMBER DIMENSIONS VOLUME ft. 
Slabs 

peccvom floor 1 Moree eye ew xX ore Owe 
Madale floor 1 ioe. cl esc eX ea 20,379.36 
Roof 1 Woe “2uxanic et x 6 62605205 
Supports 

incerior 32 ey Sl 1S SS ae ee. 1,136.90 
Beecverior 20 omen! Kea ot x 2a tt 906.58 
Joiners 48 Cex aia x co! | 5! ee 2s Oye) 
Total Volume DO, Seiag5 


Thus the mass is: 
Mass = M = 56,527.93 ft x 150 lb./ft-> = 8,479, 189.5 lb. 


Using the tables in the American Society of Heating, 
Refrigeration, and Air-Conditioning Engineers (ASHRAE) Guide 
and Data Book [4], to obtain the specific heat capacity of 
eomemervemGepe = 0. 156 BTU/1b.°F), a constant representing 


the building energy storage capacity was computed: 


M x Cpe = 8,479,189 lb. x 0.156 BTU/1b.°F = 1,322,753 BTU/°F 


WA 








Weonvensore nent LOSS AND GAIN <CONDUCTION) 

Conduction of energy through the building envelope is 
the largest loss term in the energy balance equation. The 
bullding envelope is defined as "All external surfaces which 
are subject to climatic impact; for example, walls, windows, 
Boor, etloor, ete" [5]. Architects and engineers have 
developed a method for determining the amount of energy 
transferred through the envelope by simplifying the heat 


flow equation down to: 


Q — fe Oe ai 


he 
it 


The total heat flow through the envelope 
‘lagged oul Gigi o Y gar 
U = The heat conduction coefficient for the 


2p 


Struetrure inesiU/nr. ft. 

DP 

RR =sine sW@rfiace area of the structure in Ft. . 
AT = The temperature difference between internal 


: SES 
and external air temperature in F. 


This formula and the method from which it is derived is 
explained in detail in the ASHRAE Guide and Data Book. A 
brief explanation of the formula will be discussed, but it 
1s by no means comprehensive. 

Around the year 1928, a theory of how heat conduction 
through a substance could be modeled was presented. This 


method was based on the analogy of a material's heat 
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weeadierer Characteristics to an electrical circuit. 


theory spurred research into determining 


validity of the theory, and in 1942, the 


Heating and Ventilation Engineers (later 


adopted a standard for testing materials. 


The 

the value and 
American Society of 
to become ASHRAE) 


The theory was 


Simple enough. The ability of a material cr substance to 


resist the flow of heat through it would be analogous to tne 


resistance in an electrical component to the flow of 


current. The temperature difference across the substance 


woud pe the driving potential, and the “current” through 


the substance would be the heat flow. The early theory was 
based on steady-state conditions, but later analysis proved 
the theory was more than adequate to model the heat flow 
under dynamic or changing conditions. The theory stated 
that if a collection of materials or substances were stacked 
together, then the total resistance to heat flow of the 
stack would be the sum of the individual resistances. [In 
the electrical analogy this would be equivalent to the 
resistances being in series, and the total resistance of an 


n-element series circuit is: 


i Sgt 


2d 


lees! 


R1 


The total resistance of the series circuit. 


ye 
ii 


The individual resistance of each element or 


ye) 
i 


material. 
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The total thermal resistance (Ry) iSude weeemoroca ll. OF 
the thermal conductance (U) and the value of (U) is the most 
important quantity to the engineer or architect in the heat 
conduction equation and in the study of an energy efficient 
building. The smaller the value of thermal conductance, the 
less heat or energy will be lost by the building for a given 
temperature difference between the inside and outside 
environment. It should be noted that this also applies when 
Ome 1S trying to cool the inside of a Structure. The calcu- 
lation of the envelope heat loss requires the calculation of 
two constants, the thermal conductance and the surface area. 
The calculation of AT requires knowledge of the internal and 
external temperature, both variables, and not precomputable. 

A single number representing the total surface area of 
moemwoulmlaing LS not useful, since for different structural 
members the thermal conductance (U) differs and each partial 
Surface area must be individually computed. Bullard Hall is 
divided into six surfaces; the roof, floor, and the four 
Walls. None of these surfaces are homogeneous, so their 
individual areas must be calculated. The following table is 


a computation for those different areas: 
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Surface 


Floor 


Roof 
Vents 
Skylights 
seuttle 


Roof area 


Walls 
Glass 
Wood 


Concrete 


Number 


Material 


concrete 


metal 
Deas cic 
metal 


tar + rock 


1/4" glass 


io one 


Dimensions 


Pee 2. iexk Sawa tt 


BAR x3 Neo 
3'5,5" X 315.5" 
3! X o'6" 


Area Pte 


ewe 0.2 


Bene 

107.64 
[se 
soe 0 


8,280.57 
3,415.64 


oe ee 


Tne thermal conductance for each type of surface must 


be calculated. 


In Appendix A a set of drawings of the 


building's surface materials and how they are assembled is 


presented. 


the thermal resistance of each structure. 


Along with the drawings are the computations for 


taote is a colPection of that information: 


iG 


The following 





Surface Thermal 


ELoor 


Roof 
Vents 

Sky lights 
Soutit lie 


Roof area 


Walls 
Glass 
wood 


Concrete 


Conductance (BTU/hr. Pt .° OF) 


0.48 


S108 
es. 
705 
St 


The only remaining computation is the complete building 


neat loss coefficient (U x A). 


The total envelope neat loss 


constant is the summation of all the areas times the 


respective thermal conductance 
the previous tables, the follo 


eaae cOomstant. 


Usameg the informat von. from 


wing table is a computation of 


1/7 





Surface Area U U x A 


Floor 17,468.02 0.48 Sse 05 
Roof 

Vents ype Bo. crs SL al 
pky lights 107.64 Te15 3) eres 
Scuttles fan 116% Pan 22 
Roof area Piese.O. 10 Os 2,946.11 
Walls 

Glass 8,280.57 Pais 9,357.04 
Wood SS., 6 4 1200 3,415.64 
Concrete 2,049.18 Oe la] 1,249.99 
Total Gaqu7nr. > 17,149.91% 


Pe liais tisure does not include the constant for heat loss 
through the floor because the floor is not subject to the 
exterior air temperature, but the ground temperature. 


B. INFILTRATION (CONVECTION) 

Convection losses in a building are caused by the flow 
of outside air through the building being heated and then 
exhausted back to the environment. This loss is called 
Mr teatron ane the cause for this loss is the pressure 
difference across the building's surfaces. These pressure 


differences are caused by two different conditions: the 
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temperature difference between the inside and outside and 
tme flow of wind across the surfaces of tthe building. Both 
of these effects cause the air to flow in and out of any 
openings in the surface of the building. A large body of 
data from the 1940's exists on these heat flow mechanisms 
but no practical model has been yet developed. This is made 
evident by the request for proposals made in 1980 by the 
Department of Energy for a suitable model of the air flow 
through the openings in a structure. A short explanation of 
these effects follows. 

When a temperature difference exists between any two 
areas there is a natural flow of air from the colder area to 
the warmer area due to the density difference of the air. 
This density difference causes the warmer air to rise and 
the colder air to settle in the lower areas. This differ- 
ence causes what is called the chimney effect. When the air 
in a building is warmer than the outside air a high pressure 
area is set up in the higher levels of the building and the 
warm air flows out of the upper levels. At the same time, a 
low pressure area is established near the floor or lower 
levels; this low pressure draws air into the building from 
ee Outside . 

lemodmcilatbes the sdifferential pressure across a building 
the engineer must first determine the building's neutral 


zone or the level above the lowest level or floor where tne 


1s, 





temperature 1S average. Once this has been determined the 
engineer can then compute the total differential pressure 
across the building. Next, the engineer must locate all 
openings in the building envelope, calculate their distance 
from the neutral zone, and compute the air flow through that 
opening because of the pressure differential. 

As air flows over a surface, high and low pressure areas 
are created whose pressure magnitudes depend on the angle 
between the wind and the surface of the structure as well as 
the velocity of the wind across the surface. This effect is 
Similar to the lift created on an airplane wing and is 
called the Venturi principle. On the windward side of the 
structure a high pressure area develops and on the leeward 
Side a low pressure area is developed. This differential 
pressure causes the flow of air through the building which 
is independent of the temperature difference that exists in 
the building. Again the theory is simple but the actual 
calculation and determination of the overall effect is 
difficult. Wind speed over the building surfaces is not 
always uniform and therefore an averaged estimate is used. 
If the exact weather conditions for a given period were 
Known, then with a reasonable amount of calculation, the 
engineer may have some idea of what actually occurred. But 
without a complete set of data, the engineer would only be 


guessing at the actual differential pressures involved. 
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An engineer or architect would not actually calculate 
Maewcirects Of the wind and the thermal difference; in fact, 
ne would employ another estimate to approximate the effects 
of both conditions. The ASHRAE Guide gives a method for 
estimating the overall effects on the building by expressing 
the neat loss in terms of the total number of times the 
entire air massS in the structure is replaced or turned over. 
For most structures, the turnover rate is one to two and a 
half times the total volume each hour, with allowances made 
Honeune Orlentation and location of openings in the 
building. Again most of the work on this method of heat 
transfer was done around the 1940's and 50's. The basic 


equation for determining the energy loss is: 


Q=V x coeff. x Cp, x p x Gye - T ) 


iLjelie 


fe 
ii 


Enecevetosso 27 BlU/hr. 
V = The volume of the structure in ft. 
Cp. = The specifie heat capacity of air 0.24 BTU/1b.-°F 


a 
ppeefedenetty. of air 0.075 1b./ft.- 


At = The external temperature in 5 
ae The internal temperature in ae 
Geese. = Ine number of turnovers in one hour. 


Since Bullard Hall is not surrounded by windbreaks or 
major wind deflecting structures, and has windows on all 


sides, the turnover coefficient for Bullard Hall is selected 
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to be ¢.5. The infiltration heat loss equation then simpli- 


mes. tO: 


Q = 454,168.22 * 2.5 * 1.08 # on = ) BTU/hr. 


Tint 
+ This heat loss or gain term is included in the model 
in ~ Sout? 


C. SOLAR GAIN (RADIATION) 

Solar power or usable energy from solar radiation is a 
new energy field. Daniel Berhrman, in his 1976 book "Solar 
Energy The Awakening Science", concludes that man is just 
now realizing the potential for solar energy. But solar 
power has been used for hundreds of years, as far back as 
the middle ages. The American Pueblo Indians were building 
Structures out of adobe clay using solar energy as their 
major source of heating, and cooling. Around the turn of 
the century, E. G. Morse studied the large mass heat storage 
system for solar energy, but not until 1976, when a 
Frenchman named Thrombe evaluated the large mass structure 
for storage and transmission of its heating properties, was 
this type of storage medium taken seriously. From Thrombe's 
studies, the theory and equations necessary for constructing 
a Thrombe Wall, which could be used to keep a structure cool 
during the daylight hours and provide heating through the 
night from stored energy were presented. The curious thing 


is that if a Tnrombe Wall were constructed in the Arizona 
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desert made of adobe clay which would provide the optimum 
amount of heating and cooling, the thickness will be nearly 
the same thickness of the abode walls constructed by the 
Pueblo Indians more than 1000 years ago. The use of solar 
energy to dehydrate fruits and vegetables for storage and 
later consumption by the Indians was noted by the first 
eoOlonists to settle in the new world, so the use of solar 
power is not new. In 1955, the International Solar 
Engineering Society was formed and held its first Solar 
Poems ysoy¥mMposium in Phoenix, Arizona. This group of 
engineers and educators preached the conservation of non- 
renewable energy Ae ee and the use of renewable energy in 
the forms of solar and wind power. Until recently no one 
Pec tnts group seriously. In 1955, Daniels and Duffie 
published "Solar Energy", a book which predicted the energy 
problems of the 1970's and recommended the study and use of 
sOlar power as an alternative to petroleum products. The 
book contains examples of solar powered water heaters and 
concentrating solar collectors for providing high tempera- 
ture water and steam for energy production. 

Our sum is a Star and a small star at that, but even so 
it produces about 3.8 x 10(23) KW of energy each hour of 
this energy. Only 1.7 x 10(14) kwh of energy reaches the 
earth's upper atmosphere and either through reradiation, 


absorption, or reflection only about a third of that energy 


es 





reaches the earth's surface [6]. Assuming 75% of the 
earth's surface is water and only 1% of the energy falling 
on the earth's remaining surface is utilized or converted 
into a usable energy source, this still amounts to 1.2 x 
10(17) kwh or better than 1600 times more energy used in the 
United States, in 1977 [7]. This amount of energy is 
substantial and represents a tremendous untapped resource. 

hiouerdings and other structures, solar tracking is not 
nearly as important as determining how much of the sun's 
ememey actually reaches the building's surface. To do this, 
one must Know where the sun is in relation to the building's 
surfaces. The U.S. Naval Observatory in Washington D.C. 
Pupbeshes the “Nautical Almanac" (H.O. 214) yearly. This 
publication contains tabulated data which with a few simple 
calculations, the position, azimuth and elevation, of a star 
or planet can be obdtained for celestial navigation, and one 
of those stars, is the sun. In order to understand how this 
data is available and how it is used a brief explanation of 
celestial navigation or astronomy is necessary. 

Celestial navigation is performed in a spherical coordi- 
nate system. The center of the earth, is the center of the 
coordinate system. The equator provides one reference 
coordinate and a circle drawn through the north pole and 
Greenwich, England provides the other reference coordinate. 


With two coordinates any position on the earth's surface can 
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be specified, for example Bullard Hall is 36 degrees 35 
Minutes and 43 seconds North Lat. and 121 degrees 52 minutes 
30 seconds West Long. [8]. This set of coordinates uniquely 


Paxes the position of Bullard Hall on the face of the earth. 


Longitude 
o f 
Greenwich 


England 


Equator 


0.0 Lat. 
0.0 Long. 





Figure 1 


Celestial Coordinate System 


ime is also an important quantity of interest and to 
Standardize time all celestial time is measured from the 
time in Greenwich, Greenwich Mean Time or GMT. 

A major assumption in celestial navigation is that all 
the light rays arriving from a star are parallel. This is a 


good assumption even for the sun since the maximum error for 
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the sun is 0.002 degrees or 8.7 sec. of are. The difference 
is hard to measure, and is of little importance in the 
calculation of solar gain. Under these assumptions, the 
position of a star can also be uniquely given by a set of 
latitude and longitude coordinates, but in the case of a 
celestial body these coordinates are referred to as the 
declination and GHA, Greenwich Hour Angle. 

With this short introduction one can now compute tne 
sun's position at any time. First, the GMT is computed. 
Maeess is done by noting the time zone (California is in time 
zone +38 or California is 8 hours behind Greenwich) and 
adding that to the local time expressed in 24-hour time. To 
compute GMT in California add 8 to the Local Time. This 
establishes the GMT. Next, one enters the Nautical Almanac 
with the date and GMT and reads the sun's GHA (Longitude) 
and declination (Latitude). From the GHA subtract the 
Longitude of the observer this is now referred to as the 
Local Hour Angle (LHA). Now applying Napier's rules for the 
solution of the celestial triangle for the sun's altitude 


ange azimuth [9]. 


AMEE LUDE = 


arcsine(sin LAT x sin DEC + cos LAT x cos DEC x cos LHA) 


AZPMUTH = 
arectan(sin LHA) / cos LHA x sin LAT - tan DEC x cos LAT)) 


26 








VERTICAL 
NORTH 
Ay i 
: RIM / 
/ a 
Ala 
7 & aa =) 1 
faa 
IN 
FIGURE 2 


Sun’s Relative Position 


27 





These two equations compute the altitude (the angular height 
above the norizon) and the azimuth (the angular measure from 
due North) of the celestial body. 

In order to have a dynamic model of the sun's position, 
for continuous computation of the energy received, it was 
necessary to find a suitable method to approximate the sun's 
Weaek COntinuous!y. This particular problem has already 
been solved by the U.S. Naval Observatory and in the 
publication "Almanac for Computers" which is published 
yearly. This publication uses a fifth-order polynomial to 
approximate the motion of the stars and planets and contains 
the coefficients to be loaded into a small computer to 
compute the position of a star at any time to within 6 
seconds of error. The procedure will not be discussed 
because of its length, but copies of the publication and 
instructions can be obtained from the Naval Observatory on 
request. 

If there were no atmosphere present, the eartn's surface 
would receive the equivalent of 429.2 BTU of energy per 
square foot each hour [7]. Because of the atmosphere and 
the chemicals and particles in the atmosphere, the earth's 
Surface receives only about 300 BTU per square foot per hour 
when the sun is directly overhead. The rest of the energy 
is apeerbed or dispersed into the atmosphere. This fact of 


Hagens referred to as diffusion and accounts for the blue 
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sky and the light we perceive as diffuse light in the 
atmosphere. This phenomenon of nature is discussed in most 
college level books on the subject of the atmosphere or 
light transmission. The actual process of how this energy 
is absorbed and then reradiated would require many equations 
ariepasges Of text Which is not worth the time to study in 
this thesis. The losses in the atmosphere still must be 
taken into account, but only in the simplest manner for our 
purposes. By uSing a device called a pyrheliograph located 
at a site not too remote from the actual site, measurements 
were taken to determine how much energy was incident on the 
Surace of “the building; Appendix B is a collection of that 
data and tne explanation of the approximations used. The 
approximations are based on the fact that the more atmos- 
phere the sun has to penetrate, the more energy is diffused. 
The amount of energy received versus height is a sinusoidal 
function and can be approximated in a computer as a look-up 
Eaolrem or a Continuous function of values. For the modeling 


orocess used here the continuous approximation: 
Need —eemak x sn ALT 


is used where ALT is the altitude of the sun above the 
horizon, and Imax is the maximum incident energy from the 


sun recorded on the pyrheliograph for that day. 
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Once the amount of energy from the sun per square foot 
1s tabulated another problem must be considered. While tne 
Sun traverses the sky the surface of the building which is 
fermied tO the sun's energy is changing, or, the total 
surface "seen" by the sun changes in relation to the sun. 
Since the energy of the sun on the ground has been calcu- 
lated in square feet of surface whicn is seen by the sun, 
some type of formulation or scheme is needed to take into 
account this "seen" surface. The surface normal is a 
mathematical construct which defines a vector which is 
perpendicular to the surface and haS a magnitude or length 
equal to the area of the surface. This vector now repre- 
sents the total surface and can be manipulated vectorially. 
To find the image of some vector on another vector the 
Operation of the dot product is used. If we consider the 
See so pesition as the direction of the sun vector and the 
amount of energy as the length of this vector we now have a 
means of determining how much energy is received or incident 
on any surface. The dot product of the sun vector with the 
surface normal vector represents the energy incident on the 
surface. The following figures represent the mathematical 
construct of the normal vector and the dot product in 
graphical form. 

The solar energy portion of the model has become 


extremely complicated at this point, and would appear to be 
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complete except for one thing; not all the energy incident 
on the building surface is absorbed by the surface. If the 
Surface were mirror-like almost all of the energy would be 
reflected back into the atmosphere. If the surface were 
true black, almost all the energy would be absorbed by the 
surface. And if the surface was a perfectly clear glass 
most of the energy would be neither absorbed nor reflected, 
but would be transmitted through the glass. Tnese three 
examples are the extreme cases and are easy to handle. The 
problem is that most buildings have an assortment of all 
three types of surfaces, in varying degrees. This variation 
is a real problem when one considers that the ability of a 
surface to absorb or reflect light energy may depend on the 
Orientation of the surface to the light energy. Using the 
example of the glass mirror we find the following: If the 
sun is shining on the mirrored surface most of the energy 1s 
reflected. If the sun is shining on the rear of the mirror 
or leaded surface, most of the energy is absorbed. finally, 
Peamraesmirnror surface is parallel to the light, most of the 
energy is transmitted through one edge of the glass and out 
the other. Again, most of the data collected and studied is 
from the early 1940's to the late 1950's, but the use of 
selected absorbers and reflectors in the collection of 
energy has sparked new studies in this area. The basic 


energy balance equation describing the problem from [10] is: 
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a = percentage of energy absorbed by the material 

t = percentage of the energy transmitted thru the 
material 

r = percentage of the energy reflected by the 


material 


The actual calculations of reflection coefficients were not 
done because of the lack of suitable test equipment and the 
complexity of the problem of assigning values for each 
incident angle of the sun on the various surfaces. An 
approximation for the total overall absorption coefficient 
was selected to be 0.85 mostly based on reports quoted in 
mies I9o0 book “Thermal Properties of Buildings” [11]. 

The final problem to be considered by this model was to 
miemude the effects of objects blocking the incident solar 
energy. Most structures do not stand isolated by them- 
selves, but are usually grouped together causing the problem 
of shadowing. Shadowing reduces the amount of solar energy 
which eventually reaches the structure. If one wishes to 
receive the maximum amount of solar energy the structure 
must not be surrounded by objects which will shade the 
collecting surface. This shadowing may be from buildings, 
trees, or the general terrain surrounding the structure. 


Available these days are small hand held shadow angle 


38 





protractors. These devices are actually templates which 
Gmeart the position of the sun (altitude and azimuth) for a 
peective (latitude and time of day. The architect or builder 
orients the template to due south at the site where the 
BURUCtCURE 15 to be built. Then by noting the objects in the 
field of view of the template, he can determine when the 
Structure will be shaded. Using this knowledge he can then 
recommend, using the site, or finding a better location. A 
shadow angle protractor was not used in the determination of 
the solar cutout angles for Bullard Hall. Instead, the time 
the sun first struck Bullard Hall was noted. Then, by use 
of the altitude and azimuth formulas taken from the U.S. 
Almanac for Computers the "cutout angle" was calculated to 
be 15 degrees above the observer's horizon and entered into 
the program. In the morning, after sunrise, the east and 
south faces of the building are partially shaded by trees 
and Spanagel Hall, the 15 degree solar cutout takes this 
shading into account. The west side of the building is 
shaded in the afternoon by trees and hills. The cutout for 
the west side was also calculated after observation and 


determined to be 34 degrees above the observers horizon. 


D. AIR HANDLING AND CONDITIONING 
An important item which so far has been omitted from tne 
model is the air handling system. This item could be a 


Substantial part of the building's environment and one of 
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Figure 5 


Air Handling Equipment 
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the chief resources available to the engineer in controlling 
the atmosphere in the structure. Supply and exhaust systems 
contribute to the infiltration losses by causing shifts in 
the neutral zone. They also provide a means of causing more 
air to be exchanged with the environment than the total 
effect of infiltration. The heating and air-conditioning 
systems can add or remove heat, humidify or dehumidify the 
alr, and filter the atmosphere. The architect and builder 
have the most control of these systems, and these systems 
are the ones which must be considered when the optimum 
@onurol strategy is devised. 

Bullard Hall is an extremely simple system with respect 
to air nandling systems. It contains only three exhaust 
fans and only one of these fans ig normally on, the fan that 
takes a suction on the lavatory facilities. This fan 
exhausts nine hundred cubic feet of air per minute (CFM), 
and has already been considered in the infiltration losses. 
The other two fans (see Figure 5) are main exhaust fans, 
capable of removing 19,400 CFM from the building. Although 
installed, operational, and controllable through the 
Honeywell AD-1000 heating control system, the fans have not 
been on except for routine testing since the control system 
was installed. Therefore they were not considered in the 


dynamic model. 
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Bullard Hall has no air conditioning system but does 
Seacaives/ Unit heaters dispersed about the first and second 
meQOnS.. inese unit heaters are controlled individually or 
in groups of three by 25 thermostats placed throughout the 
building. <A central thermostat controls the valve of the 
main steam supply to the building. The Honeywell system 
Operates the steam supply valve opening it at preselected 
times and temperatures. If the temperature measured by tne 
central thermostat on the lower floor of Bullard Hall 
indicates less than 65 degrees Fahrenheit and it is between 
0600 and 1600 on a weekday the steam valve will open and 
allow steam into the building up to the heaters. If the 
local heater thermostats sense a temperature below their 
setpoints, the unit heaters will come on and heat the 
building. This on/off function is easily modeled as a step 
function which is activated when the internal temperature 
drops below 65 degrees and is zero when the temperature is 
above the 3-degree dead zone set into the thermostats. fhe 
heaters in Bullard Hall are capable of delivering 1,937,600 
BTU/hr. which is equal to the maximum expected neating load 
on the building at the design temperature difference 
(internal temperature of 72 degrees and the outside 


temperature of 36 degrees). 
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Eee mhorce AND MACHINERY 

People and machinery in the building also add to the 
heat and cooling load, they usually provide such a small 
amount of energy as compared to the other sources that they 
are neglected except in unusual cases, such as high-density 
office buildings and computer operating spaces of buildings. 
The human body gives off an average of 400 BTU/hr. [4]. The 
major source of machinery heat in buildings, is in the form 
of lighting. (A 100-watt light bulb gives off about 340 
BTU/hr.). Bullard Hall is seldom occupied by more than 
fifty people and the maximum number of lights noted during 
the period 1 June to 1 December 1980 was 85. Neither of 
these sources was ever significant though they are included 
in the model. 

The biggest difficulty encountered in the modeling 
problem was with the people, but not that their presence 
could not be accounted for in the model. The problem was 
that people would leave the exterior doors on opposite sides 
of the building propped open which created a tremendous 
draft and more than doubled the infiltration losses. fTnis 
problem was so acute that the only reliable and consistent 
data taken was when the building was nearly empty and the 


entrances could be watched. 
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Itl. COMPUTER MODEL SIMULATION 


The digital model is one of the means available to the 
engineer to test possible modifications of a system without 
actually investing the resources in a modification. This is 
not to say that the model will truly simulate the system 
under all conditions, but with a good model, problems in a 
system can be spotted early and costly changes might not 
be necessary. The IBM 360-computer using the Digital 
Simulation Language package provided with the machine was 
used aS the method for simulating the building. This system 
Was chosen because the language is simple and easy to 
annotate, and the IBM 360-computer is normally the only one 
available to the student. 

ine followimg flow graph (Figure 6) is a pictorial 
representation of the computer programming used in the 
medelings of the building. The heart of the process is the 
integrator network where the output (TINT) represents the 
building's internal temperature. The driving force for a 
change in that internal temperature is lumped together in 
CRG.  inis quantity is the Q or rate of energy transferred 
to the system as heat in the energy balance equation derived 
tieedemmogel echaprer> OQOTOT also contains the infiltration 


term MCe,. - eee) derived in section II.B. of the model 
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Ceeocer. @me driving force is divided by the mass and 
specific heat of the building material and then entered 4s 
the function to be integrated. The resulting outpnut (TINT) 
is then used to recompute the temperature difference (DLT) 
between the internal and external (TEXT) air temperature and 
Gm@eweround loss Bbemperature difference. An error calcula- 
tlon is then made between the actual internal temperature 
and the calculated value; then the important values are 
Caleulated and displayed in grapnical form. A& complete copy 
Of Ome programas run for 18 October 1980 is contained in 
Meecnaix C. Ihe comparison of Ghis model with other 
ek@ocimey models, particularly the Building Load Analysis 
System Thermodynamics model, was not done because the BLAST 
model arrived too late to be installed and become opera- 
Piemal at the computer facility. It Ws hoped that work 
Might be done with the BLAST simulation program in future 
thesis studies. 

The following set of figures represents the outcome of 


Ts 
ta 


the medeling effort. Figure (7) shows the actual and 
Cabeulaved temperatures in the building. Figure (8) is the 
graph of the absolute error between tne two temperatures. 
Figure (9) is a piecewise linear graph of the actual 
external temperature surrounding the building. Figure (10) 


is a representation of the difference between the internal 


temperature (measured) and the actual outside air 
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temperature. Finally, Figure (11) is the calculated total 
Geaving force (O101) om the building from all sources. 

It can be noted that in Figure (8) that the maximum 
difference between the actual and calculated temperatures 
Was about 2.5 degrees Fahrenheit. This difference is mostly 
due to the heating effects of the sun as indirect heating 
ameenOG as Sunlight actually striking the building. This 
effect was not included in the model and no reference could 
be found that even hinted as to how this effect could be 
imeluded in the model. This particular simulation was the 
best of a series done for different days and different 
weather conditions. This particular model works best on 
Ccetewor MaZy days. The largest errors occurred on days 
when cioud cover or rain made it almost impossible to 


determine when the sun was actually on the building. 


moueecona ROL AND COST 

The object of the model was to emulate the character- 
ficmresmonetne building, so that it would be possible for the 
engineer to devise a control scheme. The control scheme 
would provide the maximum comfort for the individuals in the 
bullding, while minimizing the amount of fuel necessary to 
keep the building comfortable. In order for the engineer to 
devise the optimal control scheme he must place a price on 
eacn of the parameters he wishes to optimize. This price is 


Pererred to as the "cost" of that particular parameter, the 
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cost of fuel is easy to place a value on, but the "cost" of 
human comfort while occupying the building is hard to 
determine. In 1979, the Naval Postgraduate School spent 
over $240,000 for natural gas, which is used as the primary 
fuel for heating. In 1980, the price has more than doubled. 
Memron selecting a cost for the comfort of the individual, 
the health community has still been unable to place an 
estimate on how the performance of people is affected by the 
environment except under extreme working conditions. The 
engineer is faced with the problem of selecting a cost and 
then of selecting the best cost function to optimize. The 
Auer Lye co mprek thewcorrect costs is more of an art than a 
sclence, and 1S based mostly on intuition and experience. 
The present control theory used on the Honeywell system is 
one based on a minimum time response, when the temperature 
in the building reaches a preset temperature of 65 degrees 
the system immediately allows the maximum amount of 
Steam/energy available to flow into the building, thus 
raising the building's internal temperature back Ze) a 
desired level (07 degrees). The only modifications to that 
type of control system nave been made by overriding the 
Saesoutetunebion of the orogram during specific times of the 
day or night. This override reduces the energy used by the 
Digwidtimgmwaile it 1s not occupied, but it may not be the 


most energy efficient method to heat the building. 
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Pitseeent cCOst functions” wili produce different types 
of control functions. The minimum time response cost 
function is but one of many cost functions which might be 
eonsidered for the control of Bullard Hall's heating system. 
Another cost function to be considered could be the minimum 
fmwel cest function. This cost function attempts to meet the 
system requirements (maintaining the temperature in the 
building) while expending the least amount of fuel. Again 
insufficient time was available to develop any new control 
scheme for the temperature control in the building using the 
cost functions previously mentioned. One other factor which 
previously had not been considered in the control scheme is 
the use of the ventilation fans as a control mechanism. [f 
the fans on the roof were connected by an air-handling 
system and used to recirculate the air in the building as 
well as just exhaust the internal air, the building's 
internal environment would improve due to the cooling of hot 
Spots (the south facing wall) and heating the colder areas 
in the building (Room 112). These ideas again are only 
suggestions in the nope that someone else may attempt to use 
the existing knowledge and computer models (particularly the 
BLAST model) and pursue a more energy efficient control 


program. 
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IV. DATA COLLECTION AND INSTRUMENTATION 


The main objective here was to develop a working 
thermodynamic model of a building and verify that the model 
did indeed simulate the building. The problems of 
developing the model were far outweighed by the problems 
encountered in the collection and verification of data. The 
main parameter of interest, the internal building tempera- 
ture was measured hourly for as many as three days at a 
Cime. Since remote temperature sensing and storage devices 
were not available, the data was collected by hand and 
recorded. This process alone was very time consuming. If 
further investigation into the thermodynamic effects of the 
environment are to be conducted, it is recommended that the 
first order of business be to build a set of remote sensing 
devices and tie them together to a microprocessor so that 
the data can be collected continuously. 

The second biggest problem in the modeling process was 
to determine the amount of solar energy that was incident on 
the structure. A pyrheliograph located at the Hopkins 
Marine Station in Pacific Grove was used to estimate the 
amount of solar energy that was available at the building 
under study. The instrument was approximately four miles 


from the site and could only really be used aS a measuring 
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device when there were no clouds in the sky or when the 
total area was covered with fog. At one time the Naval 
Postgraduate School had an operational pyrheliometer but 
because of age and lack of funds the instrument has been 
allowed to become non-operational. This particular instru- 
ment would be advantageous to the Navy and the civilian 
Semmunicy to Nave functional, since so little information is 
avallable about the amount of solar energy actually avail- 
able on the peninsula. Recently the U.S. Army contracted to 
build over 300 homes at Fort Ord that will use solar panels 
aS a booster for the hot water systems in the houses. The 
Army engineers who are responsible for the evaluation admit 
Diete tune, Ouclay of funds to install the solar systems may 
never be recovered by using the solar panels, but the 
information that they hope to gain concerning the area and 
the availability of solar energy in the area is considered 
to be worth the expense. Presently the Meteorology 
Department at the school provides the local news agencies 
with most of the temperature and meteorological data on the 
area. The Fleet Numerical Weather Center, which is not part 
of the school command, contains the facilities to evaluate 
almost any set of atmospheric conditions for the area with 
their facilities, yet they remain untapped. 
Waewinstrumentation availability and the collection of 


datamarte sre most difficult tasks in developing the 
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Simulation. Because of the upcoming computer replacement 
and the large volume of computer programs that were run at 
Gees school during the last three months, little time could 
be spent in trying to develop a more complete or reliable 
model. The BLAST model is considered by many in the 
engineering field to be the most accurate and complete 
System to study the effects of different control schemes on 
a structure. Since BLAST is incompatible with the present 
computer system and required the use of the computer 
facilities at the Weather Center, little testing or 


comparison of the two models was done. 
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V. CONCLUSION 


The actual model developed is extremely simple and works 
very well under ideal conditions, even though the model 
fails to take into account the indirect solar gain caused by 
reflections from other buildings and diffused light. The 
model, though incapable of handling partially cloudy or 
rainy days, is an aid for those people who have no concep- 
tion of how the temperature inside a building will react to 
the outside environment. The model, though simple, contains 
all of the major forces acting on the building and can be 
used in developing a better control scheme for the tempera- 
tire control’ in Billard Hall. 

The major reason for using the model was to gain 
Anowledge concerning conservation and more effective use of 
solar energy. Presently, the U.S. Navy and the Naval 
Research Laboratory have more facilities and programs aimed 
at conserving petroleum products and uSing renewable energy 
than any other branch of government except the Department of 
Pitemsye. yet the officers at this institution must look 
elsewhere for information on energy programs and conserva- 
tion measures. A large portion of this school's graduates 
go on to jobs where it would be possible to implement energy 


conservation measures and save the Navy large sums of money, 
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but few have ever been exposed to the information they need. 
Information on conservation and alternative energy sources 
1s not readily available and few students have pursued any 
courses of study dealing with energy conservation. 

Most of the structures at the Naval Postgraduate School 
are not designed to be energy efficient, but with research 
efforts by students and faculty many of the buildings could 
be modified and used to demonstrate energy saving ideas and 
methods. Bullard Hall, as can be ascertained, is a 
Staucture that wasties energy. The building has too many 
doors, some of them never used. The thermal circulation 
in the building is inefficient, yet no one has suggested a 
modification to the present ventilation system which now 
exhausts the warmest air, even when shut off. Room ile 
1S a proverbial ice box; the sun never shines into the 
Space and the temperature is consistently five to ten 
degrees colder than the rest of the building. The room 
is glass on two sides and the blinds are almost always 
closed. This prevents little reflected light from entering 
the back of the room, but allows energy to escape. The 
meom also contains a set of doors that are never used. 
These doors do not seal tightly and the room is perpetually 
digaiiviae 

One issue this work did uncover was the question of 


whether sufficient instrumentation existed at the boiler 
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plant to actually measure the energy being converted in the 
boilers and where this energy goes. Few of the readings 
taken by the boiler plant operator have any meaning and much 
of the instrumentation is either out of calibration or 
inoperative. An example of this is the house steam 
differential pressure detector. In May of 1979 something 
either broke in the instrumentation circuits or the 
differential pressure cell became inoperative due to an 
Seseruccroen Ot a part failure in the cell. Either way, the 
recorded steam demand immediately dropped to one half of 
wnat it had been, without any change in any of the other 
measurable parameters. To the author's Knowledge the 
instrument is still out of service and no attempt has been 
made to check the problem. 

Suggestions for future study and investigation include: 
1) A joint Mechanical and Electrical Engineering Department 
group deSigning and constructing a microprocessor based 
remote temperature monitoring and data collection system. 
2) A Mechanical Engineering Department study and cost 
analysis of installing insulation in isolated rooms and 
energizing air-handling equipment in Bullard Hall. 3) A 
repair of the steam plant instrumentation and other energy 


monitoring equipment at the Naval Postgraduate School. 
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APPENDIX B 


PRYHELIOGRAPH DATA AND IMAX ASSUMPTIONS 


The following chart represents a complete data week from 
the Hopkins Marine Stations pyrheliograph installed in 
Pacific Grove. By using the calibration procedures in the 
Belfort Instrument Manual which accompanied the instrument a 
base line measurement of 430 BTU/ft. hr was established at 
approximately 12.5 mm. above the zero input line. Using 
this base line as a reference the maximum solar intensity 
for a particular day can be computed by forming a ratio of 
the maximum neight as compared to the 12.5 mm. reference 
mime. A pictwremor the instrument is also included for 
information. 

Using the data collected for September and October 1980 
Was compared to the expected values for the area as recorded 
in the ASHRAE Guide and Data Book, results were within 8% of 
the expected value quoted for the area. In the model the 
parameter Alpha represents the maximum incident solar energy 
Reeerded for that day times 0.85. Phe 0.85 correction 
figure was chosen to indicate the amount of energy incident 
on the building which is absorbed. This approximation is 
based on data taken from the Reference 11 "Thermal 
Properties of Buildings" on the estimated absorption 


coefficient for concrete. 
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